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ABSTRACT 

We present observations of the Rossiter-McLaughlin effect for two exoplanetary systems, revealing the 
orientations of their orbits relative to the rotation axes of their parent stars. HAT-P-4b is prograde, with 
a sky -projected spin-orbit angle of A = -4.9 ± 11.9 degrees. In contrast, HAT-P-14b is retrograde, with 
A = 189.1 ±5.1 degrees. These results conform with a previously noted pattern among the stellar hosts of 
close-in giant planets: hotter stars have a wide range of obliquities and cooler stars have low obliquities. This, 
in turn, suggests that three-body dynamics and tidal dissipation are responsible for the short-period orbits of 
many exoplanets. In addition, our data revealed a third body in the HAT-P-4 system, which could be a second 
planet or a companion star. 

Subject headings: planetary systems — planets and satellites: formation — planet-star interactions — stars: 
rotation 



1. INTRODUCTION 

The Rossiter-McLaughlin (RM) effect, a spectroscopic 
phenomenon that occurs during stellar eclipses, has recently 
been used to study spin-orbit alignment for transiting exoplan- 
ets. Although the first 9 published results suggested that the 
orbits of close-in planets are all well-aligned with the equa- 
torial planes of their parent stars (Fabrycky & Winn 2009), 
the next 20 results were more diverse, including orbits highly 
inclined with respect to the star's equatorial plane (see, e.g., 
Hebrard et al. 2008, Winn et al. 2009a, Johnson et al. 2009) 
and even retrograde orbits (Anderson et al. 2010, Narita et 
al. 2009, Winn et al. 2009b; Triaud et al. 2010). 

These results have been marshalled as evidence against 
the standard scenario for planet migration, in which disk- 
planet tidal interactions cause the planet to spiral inward. In- 
stead the results suggest that many close-in giant planets ar- 
rived at their current locations through gravitational pertur- 
bations from other massive bodies, followed by tidal dissi- 
pation (Triaud et al. 2010, Winn et al. 2010a, Matsumura 
et al. 2010). Another possibility is that protoplanetary disks 
are frequently misaligned with the rotation of their host stars 
(Bate et al. 2010, Lai et al. 2010). 

Recently, a possible trend emerged from the results: mis- 
aligned systems tend to have stars with effective temperatures 
exceeding about 6250 K, or masses ^1.2 M Q . The evidence 
for this pattern is based not only on RM observations (Winn 
et al. 2010a) but also on the line-of-sight stellar rotation ve- 
locities of transit hosts (Schlaufman 2010). This trend may 
indicate that planet formation and migration are fundamen- 

1 Department of Physics, and Kavli Institute for Astrophysics and Space 
Research, Massachusetts Institute of Technology, Cambridge, MA 02139 

2 Department of Astronomy, University of California, Mail Code 341 1, 
Berkeley, CA 94720 

3 Department of Astrophysics, and NASA Exoplanet Science Institute, 
California Institute of Technology, MC 249-17, Pasadena, CA 91125 

4 Las Cumbres Observatory Global Telescope Network, 6740 Cortona 
Drive, Suite 102, Santa Barbara, CA 93 1 17 

5 Department of Physics, Broida Hall, University of California, Santa 
Barbara, CA 93106 

6 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cam- 
bridge, MA 02138 



tally different for low-mass stars than for high-mass stars, for 
which there is already evidence in the distributions of planet 
mass and period (Bowler et al. 2010). Another possibility is 
that the formation and migration processes are similar, but that 
the subsequent tidal evolution is different (Winn et al. 2010a). 
In this hypothesis, cool stars are observed to have low obliqui- 
ties because tidal evolution drove them into alignment, while 
hot stars retain their "primordial" obliquities because of their 
thinner (or absent) outer convection zones and consequently 
slower rates of tidal dissipation. 

Although the trend seems clear, it is difficult to assess its 
true significance because many possible variables were ex- 
amined before alighting on stellar temperature and mass. The 
only way to be sure is to gather more data. This paper presents 
results for the next two systems we observed after the trend 
had been identified. Both systems have short-period giant 
planets, but HAT-P-4 is "cool" (T e ff = 5860 ± 80 K; Kovacs 
et al. 2007) while HAT-P-14 is "hot" (6600 ± 90 K; Torres et 
al. 2010). We present the observations of these systems in § 2, 
the analysis and results in § 3, and a discussion in § 4. 

2. observations 

Our spectroscopic observations employed the High Reso- 
lution Spectrograph (HIRES; Vogt et al. 1994) of the Keck I 
10m telescope, on Mauna Kea, Hawaii. We gathered 35 spec- 
tra of HAT-P-4 on the night of 2010 March 29/30, and 44 
spectra of HAT-P-14 on the night of 2010 April 27/28, in both 
cases spanning a predicted transit of the planet. An additional 
14 spectra of HAT-P-4 were gathered on other nights, at es- 
sentially random orbital phases. 

We used the standard instrument settings and observ- 
ing procedures of the California Planet Search (Howard et 
al. 2009). The iodine gas absorption cell was used to track 
the instrumental response and wavelength scale. The relative 
radial velocity (RV) of each spectrum was measured with re- 
spect to an iodine-free template spectrum, using a descendant 
of the algorithm of Butler et al. (2006). Measurement errors 
were estimated from the scatter among the fits to individual 
spectral segments spanning a few Angstroms. Tables Q] and [2] 
give the RVs, including re-reductions of data presented earlier 
by Kovacs et al. (2007) and Torres et al. (2010). 
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We also conducted photometric observations of HAT-P- 
4, in order to refine the transit ephemeris and other system 
parameters. We observed the transit of 2010 March 29/30 
(the same night as the Keck observations) with the Faulkes 
Telescope North (FTN) 2m telescope, on Mauna Haleakala, 
Hawaii. We used the Spectral Instruments camera with an 
SDSS i filter and 2x2 binning, giving a pixel scale of 0"304 
and a 10'5 field of view. Unfortunately the guider malfunc- 
tioned that night. The transit of 2010 May 7/8 was observed 
with the Fred L. Whipple Observatory (FLWO) 1.2m tele- 
scope on Mt. Hopkins, Arizona. We used KeplerCam with an 
SDSS i filter and 2x2 binning, giving a pixel scale of 0"67 
and a 23' 1 field of view. After standard debiasing and flat- 
fielding operations, we performed differential aperture pho- 
tometry of HAT-P-4 and several other stars in the field. 



3. ANALYSIS 

3.1. HAT-P-14 

We begin with HAT-P-14, for which the analysis proved 
simpler. The model for the RV data took the form 



Vcalc(0 = ^oi-b(f) + V RM (0 + 7, 



(1) 



where V or b(?) is the radial component of a Keplerian orbit, 
VrmW is the anomalous velocity due to the RM effect, and 
7 is an arbitrary offset related to the barycentric RV of the 
template spectrum. To model the RM effect, we used the 
technique of Winn et al. (2005), which entails the construc- 
tion and Doppler analysis of simulated spectra exhibiting the 
RM effect. The resulting formula for the anomalous velocity 



V RM (t) = A/(f) Vp (f) 



1.58-0.883 
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8.4kms 
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where A/ is the fractional loss of light during the transit, v p 
is subplanet velocity (defined as the line-of-sight component 
of the stellar rotation velocity at the position of the photo- 
sphere directly behind the center of the planet), and the fig- 
ure of 8.4 km s _I is the value of vsin/* estimated by Torres 
et al. (2010). In calculating A/, we adopted a linear limb- 
darkening law with a coefficient of 0.6288, based on interpo- 
lation of the tables of Claret (2004). In calculating v p (t), we 
neglected differential rotation, and allowed the stellar rotation 
axis and the orbit normal to be separated by an angle A on the 
sky plane. For a diagram of the coordinate system, see Ohta 
et al. (2005) or Fabrycky & Winn (2009). 

Many of the parameters of the Keplerian orbit and of the 
eclipses have been tightly constrained by previous obser- 
vations. We adopted priors on those parameters, to keep 
the number of free parameters in our model to a minimum. 
Specifically, we used a fitting statistic 
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where the first term is the usual sum of squared residuals, and 
the other terms enforce Gaussian priors. In this expression, 
^ciays is the orbital period in days, T c is a particular time of 
inferior conjunction (in the HJDutc system); 7d ay s is the time 
between first and fourth contact; T d a ys is the time between first 
and second contact; Rp/R* is the planet-to-star radius ratio; 
K+ is the radial velocity semiamplitude of the star's Keplerian 
orbit; and vsinz* is the line-of-sight component of the star's 
equatorial rotation velocity. All the numerical values and un- 
certainties are taken from Torres et al. (2010). We held con- 
stant the orbital eccentricity e = 0.107 and argument of peri- 
center oj = 94°, since those parameters have little effect on the 
model once the transit ephemeris is specified. 

For oy, we used the quadrature sum of the internally- 
estimated measurement error (typically 4 ms" 1 ; see Table 1) 
and a "jitter" term of 7.8 m s" 1 , a value determined from the 
condition \ 2 = AW- Thus, the excess noise was assumed to 
be uncorrected in time. It is comparable in magnitude to the 
jitter term of 7.3 m s" 1 used by Torres et al. (2010). 

There were only two completely free parameters in our 
model: 7, the overall RV offset; and A, the projected spin- 
orbit angle. Parameter optimization and error estimation were 
achieved with a Markov Chain Monte Carlo (MCMC) algo- 
rithm, using Gibbs sampling and Metropolis-Hastings step- 
ping. Table[3]summarizes the results. The first section lists all 
of the adjustable model parameters, for which uniform priors 
were adopted. The last section gives some results for other 
quantities that were computed based on the model parame- 
ters, or taken from Torres et al. (2010). The quoted values and 
ranges are based on the 50%, 15.85%, and 84.15% confidence 
levels of the marginalized posteriors. Figure Q] shows the RV 
data and the results for vsin/* and A. 

The result for A is 189.1 ±5.1 degrees, indicating that the 
directions of orbital motion and stellar rotation are nearly op- 
posite as projected on the sky. This result could be anticipated 
from a visual inspection of Figure Q] which shows that the 
anomalous RV was a redshift in the second half of the transit 
and (less obviously) a blueshift in the first half. This is an 
inversion of the more familiar pattern of a well-aligned sys- 
tem. The orbit of HAT-P-14b is strongly misaligned with the 
rotational plane of its parent star. 

Although the finding of a retrograde orbit is robust, the 
small uncertainty in A depends critically on the prior con- 
straint on vsin/*. If that constraint is dropped, a slower ro- 
tation rate is favored (vsim* = 4.5 ±2.4 km s" 1 ), allowing a 
broader range of spin-orbit angles (A = 192. O^ 6 degrees). 
These results are also illustrated by the dotted lines in Fig- 
ure Q] One possible reason for the smaller result for vsin/* 
is differential rotation: the RM effect depends on the rotation 
rate over the range of latitudes spanned by the transit chord, 
which may differ from the spectroscopically-estimated equa- 
torial rotation rate. The near-grazing transit of HAT-P-14, in 
particular, may produce the most extreme possible difference. 
Further observations with a higher signal-to-noise ratio might 
be able to identify the specific signal of differential rotation 
(Gaudi & Winn 2007). 
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FIG. 1. — Results for HAT-P-14. Left. — Apparent radial velocity variation on the night of 2010 April 27/28, spanning a transit. The top panel shows the 
observed RVs. For the bottom panel, the best-fitting orbital model was subtracted, thereby isolating the anomalous RV due to the RM effect. The black curve 
shows the best-fitting model with a prior constraint vsini* = 8.4 ± 0.5 km s _1 , and the dotted curve is the best-fitting model with no prior constraint on vsini*, 
Right. — Joint constraints on A and vsini*. The contours represent 68.3% and 95.4% confidence limits. The marginalized posterior probability distributions are 
shown on the sides of the contour plot. The solid and dotted curves show the results with and without the prior constraint on vsini*. 



3.2. HAT-P-4 

The case of HAT-P-4 was more complicated, partly because 
our RV observations revealed a third body in the system. Fig- 
ure |2] illustrates that a single planet on a circular orbit no 
longer provides a satisfactory description of the data. The 
best-fitting model has x 2 = 925.7 with 18 degrees of freedom. 
The residuals have an rms of 15.9 m s" 1 and are highly corre- 
lated, with almost all of the most recent RVs lying above the 
model curve. Allowing the orbit to be eccentric reduces x 2 to 
906.0, but the residuals still have an rms of 15 m s" 1 and show 
the same pattern. 

A much better fit is obtained when the circular orbit is 
supplemented by a constant acceleration 7. This model, il- 
lustrated in the lower panel of Figure [2] gives \ 2 = 137.0 
and more randomly scattered residuals with an rms resid- 
ual of 6.8 m s" 1 . An MCMC analysis gives 7 = 0.0246 ± 
0.0026 ms" 1 day" 1 . 

The constant acceleration may represent the RV variation 
due to a companion star or planet whose orbit is longer 
than the three-year span of our observations. Assuming it 
is a low-mass body on a nearly circular orbit, we may set 
7 ~ GM C sin ic/cfi., giving an order-of-magnitude constraint 



M c sin i c 



Jup 



10 AU 



'5.03 ±0.53. 



(4) 



However, given the limited time sampling of our data it is also 
possible that the companion has a shorter period. We have not 
found any compelling two-planet models, but we are contin- 
uing to gather additional RV data and will report elsewhere 



on the results. The pertinent conclusions for this study are (1) 
there is a third body in the system, and (2) the mass and orbital 
parameters of HAT-P-4b are subject to systematic errors due 
to the unknown influence of the third body on the RV data. 
For the latter reason, in our analysis of the transit data we did 
not employ a prior constraint on K+, as we did for HAT-P-14. 

Another difference in our analysis is that we fitted the avail- 
able photometric data along with the transit-night RV data, 
because the new data offer better constraints on the transit 
ephemeris, depth, duration, and partial duration. Our photo- 
metric model was taken from the Transit Light Curve project 
(see, e.g., Holman et al. 2006, Winn et al. 2009c). In brief, 
we used the Mandel & Agol (2002) formulas for a quadratic 
limb-darkening law, as implemented by Pal (2008). The lin- 
ear coefficient was allowed to vary freely, and the quadratic 
coefficient was held fixed at the Claret (2004) value (0.3418 
for /-band, and 0.3395 for z-band). The out-of-transit mag- 
nitude was allowed to be a linear function of airmass, to ac- 
count for color-dependent differential extinction. The errors 
for each light curve were set equal to j3a\ , where <j\ is the 
root-mean-squared (rms) residual, and f3 accounts for time 
correlations, using the method of Pont et al. (2006) as im- 
plemented by Winn et al. (2009c). Averaging times of 10-30 
minutes were used to compute /3, giving results (in chronolog- 
ical order) of 1.01, 1.73, 1.69, and 1.01. Parameter estimation 
was performed by the MCMC method. The photometric data 
and the best-fitting models are plotted in Figure [3] 

The formula for the anomalous velocity was 



Vrm(0 = A/(0 v p (t) 
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Fig. 2. — Evidence for a third body in the HAT-P-4 system. Top. — 
Relative RV as a function of time. Squares represent data from 2007, and 
circles represent data from 2010. Below the data are the residuals between 
the data and the best-fitting model involving a single planet on a circular orbit. 
Middle. — Same, but plotted as a function of the orbital phase of the planet, 
in days. The residuals from 2010 are systematically higher than those from 
2007, which is evidence for an excess radial acceleration and hence an addi- 
tional gravitating body in the HAT-P-4 system. Bottom. — Same as the middle 
panel, but for a model which also includes a free parameter 7 representing a 
constant acceleration. The result was 7 = 0.0246 ± 0.0026 m s _1 day" 1 . 
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FIG. 3 . — Transit light curves of HAT-P-4. (a) and (b): FLWO z-band data, 
from Kovacs et al. (2007). (c): FTN /-band data, from our new observations. 
These data were obtained simultaneously with the Keck RM observations, 
(d): FLWO /-band data, from our new observations. The best-fitting model 
light curves are overplotted. 



and the fitting statistic was 
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using terminology similar to that of Eqn. (0. The only free 
parameters that were wholly dependent on the RV data were 
7, A, and K+. The jitter term was ay = 5.2 m s" 1 . 

Table |4] gives the results, based on the 15.85%, 50%, and 
84. 15% confidence levels of the marginalized a posteriori dis- 
tributions. Figure @] shows the RV data, and the posteriors 
for the RM parameters vsin/^ and A. In particular, the pro- 
jected spin-orbit angle is A = -4.9 ±11.9 degrees, consistent 
with good alignment between the rotational and orbital angu- 
lar momentum. 

As before, the quantitative results hinge on the prior con- 
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FIG. 4. — Results for HAT-P-4. Left. — Observed radial velocity variation on the night of 2010 March 29/30, spanning a transit. The top panel shows the 
observed RVs. For the bottom panel, the best-fitting orbital model was subtracted, thereby isolating the anomalous RV due to the RM effect. The solid black 
curve shows the best-fitting model with a prior constraint vsini* = 5.5 ±0.5 km s . The best-fitting model with no prior constraint on vsin;'* is plotted with 
a dotted curve, although it is hard to distinguish from the solid curve. Right. — Joint constraints on A and vsini*. The contours represent 68.3% and 95.4% 
confidence limits. The marginalized posterior probability distributions are shown on the sides of the contour plot. The solid and dotted curves show the results 
with and without the prior constraint on v sin U ■ 



straint on vsin/*. By repeating the analysis with no such con- 
straint, we find vsin/* = 6.4^'] km s" 1 and A = -1 degrees. 
These much broader results are also illustrated by the gray 
lines in Figure |4] The only well-constrained combination of 
those two parameters is vsin/^cos A = 5.77 ±0.41 km s" 1 . 



4. DISCUSSION 

It might seem surprising that tighter bounds on A were ob- 
tained for HAT-P-14 than for HAT-P-4, given that the signal- 
to-noise ratio of the RM effect is higher for HAT-P-4. This is 
a consequence of the difference in the impact parameter (the 
minimum sky-projected distance between the planet and the 
star, in units of the stellar radius). The interpretation of the 
RM signal is most robust for systems with a high impact pa- 
rameter, because in such cases the two key parameters vsini* 
and A play distinct roles: vsin/^ controls the amplitude of the 
signal and A controls its shape (i.e. the phase of the transit 
when the anomalous RV switches from positive to negative, 
or vice versa). For systems with low impact parameters, such 
as HAT-P-4, the shape of the RM signal is nearly independent 
of A and both parameters control the amplitude. This leads to 
a strong degeneracy between those two parameters (Gaudi & 
Winn 2007). 

Therefore, an external constraint on vsiru* is crucial for the 
determination of A in systems with low impact parameters. 
In our study we have used a prior based on the line broad- 
ening observed in the star's optical spectrum. However, it 
must be acknowledged that the resulting estimate of vsin^ is 
subject to systematic error due to uncertainties in the compet- 
ing effects of macroturbulence and other broadening mecha- 
nisms. This is especially problematic for cool, low-mass stars 



for which turbulent and instrumental broadening exceed rota- 
tional broadening; examples are TrES-1 (Narita et al. 2007) 
and WASP-4 (Triaud et al. 2010). For HAT-P-4 the situation 
is better because rotational broadening is expected to be at 
least as important as turbulent broadeningQ 

We turn now to the question posed at the beginning of this 
paper: do hot stars have high obliquities? Specifically, among 
the host stars of close-in planets, are those with T e s ^ 6250 K 
more likely to be aligned with the planetary orbits than hotter 
stars? We have found that HAT-P-4 is a well-aligned cool 
star (A = -4.9 ±11.9 deg, T ef[ = 5860 ± 80 K), and HAT-P- 
14 is a misaligned hot star(A= 189.1 ±5.1 deg, r eff = 6600± 
90 K). Therefore, these new data strengthen the trend that was 
observed by Winn et al. (2010a). 

Schlaufman (2010) found a similar pattern using a different 
technique, involving a comparison between the observed and 
expected line-of-sight rotational velocities of the stars with 
transiting planets. He described the pattern in terms of stellar 
mass rather than effective temperature. Indeed, those two pa- 
rameters are strongly correlated for dwarf stars, with scatter 
due to metallicity and age. For HAT-P-4 and HAT-P-14, the 
stellar masses are I.X^^Mq and 1 -386 ± 0.045 M© respec- 
tively. For these systems there is a clearer contrast in effective 
temperature (6.1a) than mass (1.7cr). The very different re- 
sults for A suggest that effective temperature is more closely 
related to obliquity than mass. This in turn would support the 
hypothesis of Winn et al. (2010a) that the differing obliquities 
are a consequence of differing internal structure of the stars, 
and specifically the depth of the outer convective zone, since 

7 Given HAT-P-4 's approximate spectral type of F7/G0, the expected 
macroturbulent velocity is «4.5 m s (Gray 2008, p. 443) as compared to 



the inferred vsini* of 5.5 ±0.5 m s 
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this structural difference is more closely related to effective 
temperature than mass (Pinsonneault et al. 2001). 

It will be interesting to examine all the systems for which 
the RM effect has been measured, to see whether temperature 
or mass is more important, and whether there are other vari- 
ables related to obliquity. We defer such a study for the fu- 
ture, to allow the sample size to grow substantially since the 
last such analysis by Winn et al. (2010a). One variable that 
will be especially interesting to assess is the presence or ab- 
sence of a third body. A migration mechanism involving the 
Kozai effect, which has been invoked to explain high obliqui- 
ties, requires the existence of a third body (see, e.g., Fabrycky 
& Tremaine 2007). In planet-planet scattering scenarios, a 
third body may also be present, although it could have been 
ejected (see, e.g., Chatterjee et al. 2008). Evidence for ad- 
ditional bodies has been found in some well-aligned systems 
such as HAT-P-13 (Bakos et al. 2010, Winn et al. 2010b) and 
HAT-P-4 (this study), along with some misaligned systems 
such as HAT-P-7 (Winn et al. 2009a), HD 80606 (Naef et 
al. 2001), HAT-P-11 (Bakos et al. 2010, Winn et al. 2010c), 
and WASP-8 (Queloz et al. 2010). A systematic multiplicity 
study would be illuminating (see, e.g., Narita et al. 2010), as 
would a comparison between the observed obliquity distribu- 
tion and that predicted by the Kozai model (Fabrycky & Winn 
2009; Morton & Johnson 2010). 

It will also be interesting to extend these studies to multi- 
transiting systems such as Kepler-9 (Holman et al. 2010). 
This will allow the mutual inclinations of the orbits to be de- 
termined, along with the obliquity of the star (Fabrycky 2009, 
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TABLE 1 

Relative Radial Velocity Measurements of HAT-P-14 



HJD UTC RV [m s _1 ] Error [m s" 1 ] 



2454602.85804 -187.65 3.64 

2454603.10267 -198.11 3.64 

2454603.86302 -142.90 3.70 

2454604.09555 -81.34 4.05 

2454633.99342 212.00 4.03 

2454634.93451 -93.25 3.85 



NOTE. — The RV was measured relative to an arbitrary template spectrum; only the differences are significant. The uncertainty given in Column 3 is the internal error only and does 
not account for any possible "stellar jitter." (We intend for this table to be available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance 
regarding its form and content.) 



TABLE 2 

Relative Radial Velocity Measurements of HAT-P-4 



HJDutc RV[ms 4 ] Error [ms~'] 

2454186.98523 55.67 2.44 

2454187.11242 52.15 2.12 

2454188.01161 -68.10 2.15 

2454188.07151 -75.14 1.92 

2454189.00175 -62.98 2.36 

2454189.08264 -46.17 2.18 

NOTE. — The RV was measured relative to an arbitrary template spectrum; only the differences are significant. The uncertainty given in Column 3 is the internal error only and does 
not account for any possible "stellar jitter." (We intend for this table to be available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance 
regarding its form and content.) 



TABLE 3 
Parameters for HAT-P-14 



Parameter Value 

Model parameters 

Projected spin-orbit angle, A [deg] 189.1 ±5.1 

Projected stellar rotation rate, v sin i* [kms~'] 8.18±0.49 

RV offset [ms" 1 ] 20.6 ±1.9 

Velocity semiamplitude, K* [m s _1 ] 218.9 ±5.7 
Orbital period, P [days] 4.6276690 ± 0.0000050 
Midtransit time [HJD UTC ] 2, 455, 3 14.91794 ±0.00066 

Planet-to-star radius ratio, R p /R* 0.0800 ± 0.0015 

Orbital inclination, i [deg] 83.52 ± 0.22 

Fractional stellar radius, R*/a 0. 1 127 ± 0.0033 
Other parameters (derived from model parameters, or from elsewhere) 

Transit duration, first to fourth contact [days] 0.0910 ± 0.0035 

Transit ingress or egress duration [days] 0.0294 ± 0.0029 

Transit impact parameter 0.894 ± 0.013 

Orbital eccentricity (Torres et al. 2010) 0.107 ± 0.013 

Argument of pericenter [deg] (Torres et al. 2010) 94 ± 4 

Stellar mass, M+ [M Q ] (Torres et al. 2010) 1.386 ± 0.045 

Stellar radius, i?* [R Q ] 1 .468 ± 0.042 

Planetary mass, M p [M Jup ] 2.232 ± 0.058 

Planetary radius, R p [K Jup ] 1.142 ± 0.033 
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TABLE 4 
Parameters for HAT-P-4 



Parameter Value 

Model parameters 

Projected spin-orbit angle, A [deg] —4.9 ± 1 1.9 

Projected stellar rotation rate, vsinu [km s~'] 5.83 ± 0.35 

RV offset [ms" 1 ] -5.7 ±1.3 

Velocity semiamplitude, fjms -1 ] 66.9 ± 8. 1 

Orbital period, P [days] 3.0565195 ± 0.0000025 

Midtransit time [HJD UTC ] 2, 455, 285.03216 ± 0.00073 

Planet-to-star radius ratio, Rp/R* °- 08697 ^:ooo45 

Orbital inclination, i [deg] 88.76^'j| 

Fractional stellar radius, R+/a 0.1690+°;°°^ 

Other parameters (derived from model parameters, or from elsewhere) 

Transit duration, first to fourth contact [days] 0.1775+" Mai 



Transit ingress or egress duration [days] 0.01465^[j'jj[jjSj 

Transit impact parameter 0. 1 284^ ^2 

Orbital eccentricity (assumed) 

Stellar mass, M, [M ] (Kovacs et al. 2007) 1 .26 ± 0. 10 

Stellar radius, R* [R Q ] 1.617!#jj£ 

Planetary mass, M p [M Jllp ] 0.556 ± 0.068 

Planetary radius, R„ [R Jap ] 1.367^ g|2 



